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. Abstract 

Recent studies have shown that the local environment of the Sun is a complex one with 
lO \ the presence of several supernova shock bubbles and interstellar clouds. Even within this 

environment, the Sun is moving with a velocity of about 15 km/s toward the interior of 
the galaxy. We examine the changes in the local ISM environment experienced by the Sun 
I ■ in the past and then discuss the interaction of the Sun with these different environments. 

Lastly, we enquire how this must have changed the cosmic ray environment of the earth 
since this remains one of the few methods to test the changes in the local ISM environment 
over astronomical time scales. 
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1 Introduction 

Recent studies of the local interstellar medium reveal that the Sun is in a complex envi- 
' ronment where at least three and probably more bubbles or cavities of probable supernova 

origin are interacting (see e.g. Breitschwerdt, 2001, Frisch, 1995, 1998, Maiz-Apellaniz, 
2001, Smith and Cox, 2001, and references therein). 
. ,— i \ Reviews by Frisch (1998, 1995) have mapped the local interstellar matter (LISM) from 

several observations based on back scattering and extinction of soft X-rays coming from 
nearby sources. These maps show that the local interstellar medium is largely governed 
by a collection of at least three large bubbles or super bubbles that seem to be shells of 
supernova remnants of age between 5 and 10 Million years. A large bubble in which the Sun 
seems to be sitting is referred to as local bubble. These observations have been re-enforced 
by radio observations of more recent studies based on the measurements of the turbulence in 
the nearby pulsar radio emission due to the LISM plasma (Ramesh Bhat and Gupta, 2002, 
Ramesh Bhat, Gupta and Pramesh Rao, 2001). They have modeled the inferred plasma 
distribution into a 3 component model that gives a description of the LISM derived from X- 
ray observations. Redfield and Linsky (2002) have attempted to study the UV observations 
of Fell, Mgll and Call to determine the structure of the LISM up to 100 pc. they have 
shown that combined studies of absorption features and Doppler shifts indicate a fairly non 
uniform distribution of the LISM. Further high resolution studies in ultraviolet (Shelton, 
2002) and infrared bands (Franco, 2002) from select regions show that there is evidence of 
hot but quiescent plasma in the LISM. A more detailed compilation of the H Q lines in the 
Milky Way by Haffner (2002) has shown that conditions of the warm interstellar medium 
and diffuse ionised gas suggest a rich mix of filamentary structures in the ISM. Similarly, 
first results from the FUSE mission (Moos et al., 2002) suggest that UV emission in the the 
warm component of the local bubble is homogeneous indicating an old bubble. 

Some of these shells seem to be disintegrating and forming fluffy clouds that are drift- 
ing within these shells or loops (see also Ricardo and Beckman, 2001; Seth and Linsky, 
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2001; Smoker et al., 2002). Meisel, Diego and Mathews (2002) have calculated the detailed 
trajectories of several dust grains approaching the sun and shown that they are indeed of 
interstellar origin. The dust content of the bubble have been estimated from the Ulysses by 
Ingrid and Hirshi (2002). 

These observations are in good agreement of models with the ISM that involve several 
shells in dynamic equilibrium between slowly dying supernova remnants of age of the order 
of a few million years (e.g. Smith and Cox, 2001) that are slowly disintegrating into small 
clouds. 

Donato, Manrin and Taillet (2002) have attempted to model the confinement of cosmic 
rays in these bubbles and the production of radio active nucleids in spallation inside the 
local bubble. They have shown that production rates of some radioactive nuclei such as 
10 Be, 28 Al and 36 C1 can be important markers of the spallation rates in the LISM. 

It is therefore clear that we have a fairly consistent broad picture of the local interstellar 
medium (see figure 5, Frisch 1995). Another interesting feature is that the Sun itself is not 
at rest in the local rest frame (Frisch, 1995, 1998) and has a proper motion of about 15 km/s 
toward the inner part of the galaxy. Dehnen and Binney (1998) have given a more detailed 
value of the movement vector of the sun in the local ISM using the Hipparcos data. They 
show that the best fit vector for the sun is uq = 10.00±0.36 km/s, vq = 5.25±0.62 km/s in 
the plane of the Galaxy and wq = 7.17 ± 0.38km/ s perpendicular to the plane of the Galaxy 
towards the galactic plane. However, this motion in the perpendicular plane is believed 
to be cyclic with a period of 66 million years (see also Frisch, 1998). The result is that 
the sun has been encountering a large range of ISM conditions over the last few million 
years. While attempts have been made to model the dynamics of the ISM no studies have 
been undertaken to understand the interaction of the sun under these varying conditions. 
In the next section we consider the changing environment seen by the Sun over the last 
few million years. In section 3 we consider various processes by which the sun interacts 
with the local interstellar medium and in section 4 we concentrate on the modulation of the 
average galactic cosmic ray fluxes due to these differing environments and show that the 
observed cosmic ray fluxes on earth would vary significantly due to different ambient LISM 
conditions. 

2 Position and relative movement of the Sun in the 
Galaxy 

We begin backtracking the path of the sun based on our current understanding of the 
velocity of the sun in LISM and the LISM environment. For this we refer to figure 5 of Frisch 
(1995). The sun is moving in the LISM toward the region of Galactic coordinates of (1,(3) 
= (51°, 23°)with a speed of 15.4 km s" 1 ) (Frisch, 1995)with respect to the local rest frame. 
The Scorpius-Centaurus association is a large number of massive stars. Several supernova 
explosions have occurred in this region over the last 15 million years. It is generally believed 
that the Sun is inside the supernova bubble created by supernova remnants from this region. 
Smith and Cox (2001) consider the option that the local environment is a result of two or 
three supernovae that went off in the region occupied by the Local Bubble that itself could 
have originated from the Scorpius Centaurus Association. The center of the Gould belt 
of young stars is located at a mean galactic longitude of 110° with a diameter of 200 pc 
(Poppel, 1997). The Taurus and Perseus region lies in the anti galactic center region while 
the Orion loop lies in the Galactic Longitude of I about 200°. 

We trace the motion of the sun back from this location. This is in line with the inferences 
of Bash (1986) who calculated the approximate path of the Sun over a long period. The 
details of the movement in table 1. The first column of the table gives the epoch before 
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Table 1: Solar movement in the ISM 



Time (yBP)(a) dT Region Den. (cm~ 3 ) Commcnts(b) 
Start End (y) (Temp (K)) 








1.0 


10 4 


1.0 10 4 


fluff 


0.08 
(7000) 


Size 5 pc, From SCA, 

moving in orthogonal direction to Sun 20 km/s 
B= 1.5 /iG. Highly Inhomogcneous (1 


1.0 


10 5 (a) 


3.1 


10 5 


6.0 10 3 


SN(?) 


4 (104) 


Gcminga? (2,3) 


1.0 


10 4 


1.0 


10 5 


(O.l-l)lO 5 


bubble 


0.04 (106 ) 


Expanding bubble from SCA. Age 15 My 

vel 4 10 4 km s -1 orthogonal to Sun B 1.5 /i G 

crossed the Sun 250,000to 400,000 y ago (1,3,4,5) 


1.0 


10 5 


1.0 


10 7 


(0.01-1) 10 7 


Inter-arm 


0.0005 ( 100) 


(1,6) 


1.0 


10 7 (a) 


2.5 


10 7 


3.0 10 6 


HI shock 


4 (10 4 ) 


From Gould belt (7). Duration is taken as that for 
the earlier SN.(3,6) 


2.5 


10 7 (a) 


4.0 


10 7 


5.0 10 5 


Orion belt 


10 4 ( 10) 


Size 5pc, Sun must have 5 10 5 taken yrs to cross 
at current velocity (7,8) 


4.0 


10 7 


5.0 


10 7 


(2.5-4)10 7 


inter-arm 


0.005 (10) 


(6) 


5.0 


10 7 


6.0 


10 7 


(4-5) 10 7 


Pcrsius arm 


1 (1000) 


Duration uncertain(6) 


< e 


; io 7 








(9) 



(a) For SN the start and end time are the interval between which the SN spike crossed the earth (b) Bracket 
numbers in comments refer to notes. 

References: 1 Frisch (1995); 2 Ramadurai(1993); 3 Sabalska et al(1991) 4 Egger et al(1996) 5 Van der Walt 
and Wolfendale (1988) 6 Bash (1987) 7 Clube and Nupier (1984) 



Notes to the table: 

1) SCA refers to Scorpius Centaurus association at 1 = 300-360°. About 35 SNR have been identified in 
this region (Whiteoak and Green, 1996). 

2) While Frisch (1995) claims that the evidence for Gcminga is not there, Ramadurai(1993) has claimed 
the passage of a SN shock around 30,000 yrs ago that arose from a source around the age of Geminga. 
The sun crossed the vicinity of the parent star of Geminga 10 7 yrs BP while the SN itself is 3 10 5 yrs old 
(Pavlov, 1997). The temperature and density are typical values for SNR (Harwit, 1998). The shock passage 
is calculated assuming a shock front velocity of 1000 km/s (1 pc/My) and taking an effective thickness of 1 
pc (Harwit, 1998). These parameters arc similar to those derived by Szabelska, Szabclski and Wolfendale 
(1991). Based on 10 Be data they argue for the passage of a supernova shock front near the Sun about 30,000 
yrs ago. 

3) Another SN is just arriving from the SCA association, age 11 My. 

4) The bubble is from the SN shell from the Upper Centaurus Lupus of SCO. Age 14-15 My BP. 

5) The inter-arm region is derived from figure 5 of Bertrische (1998) and Bash (1987), figure 1. 

6) van der Walt and Wolfendale (1998) argue that the gamma ray data suggest that the inter arm cosmic 
ray are of lower intensity and have a steeper spectrum. 

7) More than 130 weak line T Tauri (O type) stars have recently been cataloged in this region (for reference 
see Wichmann et al., 1997) 

8) Based on extrapolation as defined in note 5. 

9) A) We are now in the Orion structure in the inter arm region between Persius and Carina arms 

B) This is uncertain since collisions with individual objects are ignored (Bash, 1987). Clube and Napier, 

(1984) suggest that in the last 4.5 10 9 years, the sun has encountered about 10 molecular clouds. Napier 

(1985) has done a more detailed analysis to show that the Sun has had near collisions (impact parameter < 
20 pc) with 56 GMCs having > 3 10 3 M su „ and 8.2 close encounters with GMCs having M > 10 5 M sun 

C) Our galaxy is a highly evolved spiral with about 10% of its mass in diffuse gas; the rest is in stars. The 
sun is located slightly above the galactic plane, and is moving through space at 16.5 km/s with respect to 
the local standard of rest. The Sun oscillates in the galactic plane with a period of 33 My (i.e. it goes ± 56 
pc). 

D) Taking the Sun's relative velocity of 16.5 km s _1 (2 pc/My), in the local rest frame, the sun has traveled 
9 kpc in 4.5 billion years. Taking the solar distance as 8 kpc from the center, this gives a radius of 50 kpc. 
The sun therefore has traveled only about 18 percent of the distance in the local rest frame. In the galactic 
frame, the solar velocity is 250 km/s and hence the Sun has completed about 23 rotations in its life time. 
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Figure 1 : Energy density of the local interstellar medium as a function of time 



present when the Sun was in the neighborhood of a particular region and the duration for 
which it was in a specific neighborhood. It should be noted that there is a broad range of 
period when the Sun was in the inter arm region passing by the Geminga Pulsar region and 
it was inside the Orion Nebula about 25 million years ago for about half a million years. It 
may be noted that the linear interpolation of the figure 5 of Frisch (1995) would indicate the 
Sun missing the Orion arm but the arm itself has a motion orthogonal to the solar direction 
with a velocity of about 1 km/s (Nagahama et al., 1998). It seems more than a coincidence 
that it is during this period, when the local ISM was clearly more dense that the Earth 
seems to have encountered larger mass extinction periods. Since high density LISM and 
interplanetary environment would increase friction and decay teh comet trajectories quite 
considerably. However, we do not discuss this issue here. 

By extrapolation, we see that the sun passed from the Persius arm through the inter-arm 
region into the Orion belt. The sun escaped the supernova and entered the fluff in which it 
presently placed. The higher densities of diffuse gas seen within 50 pc of the sun between 
galactic longitudes 20° and 90° represent diffuse gas ablated from the parent molecular cloud 
complex. In the next section we discuss the physical processes through which the sun must 
interact with these regions before discussing the astronomical consequences. 

In figure 1 we have plotted the energy density of the local ISM as a function of time. 
On the X-axis is the time in log scale and on Y-axis we have plotted 1.5*kTisM*PiSM, 
where Tism is the ISM temperature and pism is the ISM density and k is the Boltzmann's 
constant (table 1 (see also Vahia and Lai, 2000). As can be seen from the figure, ISM 
densities have seen large changes in energy density over the last huudred million years for 
which data is given here. 



3 Physical processes of the interaction of the Sun with 
the local interstellar medium 

There are mainly three processes through which the interstellar gas is affected by stars 
(Dyson and Williams, 1997). A graphical description of these processes is given in figure 3. 

1. A massive hot star will photoionize the interstellar gas in its neighborhood and pho- 
tons of energy hv > Ih (where (/#) is the ionization potential of hydrogen) will eject 
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Stationary Ionization Front 



Figure 2: Evolution of the ionization induced interaction regions 



electrons from the atoms. The energy excess, h.v - Ih, goes directly into kinetic en- 
ergy of the detached electron. The protons are essentially unaffected by this process 
although momentum conservation demands that they experience a small recoil. The 
strong Coulomb interaction between protons and electrons produces the inverse pro- 
cess of recombination. The reversible reaction can be stated as H + hv <-> p + e~ . 
The energy of the photon produced on recombination is given by the sum of two con- 
tributions: the kinetic energy of the recombining electron and the binding energy of 
electron in the energy level into which it recombines. This binding energy is equal to 
Ih /n 2 , where n is the principal quantum number of the level into which recombination 
occurs. This is called the Stromgren Sphere. 

2. Due to photoionization, the gas temperature increases from about 10 2 K to about 10 4 
K. The ionization process itself increases the number of gas particles and therefore the 
pressure by a further factor of two. The pressure in the ionized gas is thus two hundred 
times greater than that in surrounding neutral material. This ionized gas expands and 
a shock front is set into motion. The result is the formation of an outer sphere whose 
inside edge is defined by the ionization radius and outer radius defined by pressure 
equilibrium with the local interstellar medium. This is called the ionization front. 

3. The other methods of interaction are via the solar wind (very high-speed continuous 
mass loss) which produces a mechanical piston that is bound by pressure equilibrium 
between the stellar wind and ISM. 

While all three processes occur simultaneously in the interaction of the Sun with the 
interstellar medium, in the present study we study each process independently based on the 
formulation of Dyson and Williams (1997). 

3.1 Ionization front and Stromgern sphere 

In oder to analyze the ionizing interaction, we assume that the Sun was in a homogenious 
and isotropic medium of the given environment and the transition time scales from one 
interaction region to the next were much shorter than the residence time in any specific 
medium. While this is largely a reasonable assumption, significant acceleration of cosmic 
rays and higher flux at high energies can be seen when passing through a supernova shock 
front (see e.g. Sabalska et al., (1991). We also assume that the sun had a constant luminosity 
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and produced Lyman continuum photons at a constant rate (S*). A sphere of ionized gas 
is then formed around the Sun and the boundary between the ionized and neutral material 
moves into the surroundings at a speed, which generally approaches the speed of light. The 
radius of the ionized sphere is a function of time. There are two time scales of the movement 
of the Ionization boundary or Ionization Front (IF) and the the surfaces of ionization and 
ionized gas containing the ionized gas, the Stromgrcn sphere. 

A basic property of an IF is the velocity at which it moves. If measured with respect to 
the Sun, this is a fixed frame of reference. We take the gas around the Sun be at rest (in 
the frame of reference of the Sun). At time t, the IF is at a distance R from the star and at 
time t + dt; it is at a distance R + dR (figure 3 . We take the undisturbed neutral hydrogen 
density be n m~ 3 and let J mT 2 s~ Y be the number of Lyman continuum photons falling 
normally on a unit area of IF per second. In the interior of the IF, the gas is fully ionized. 
Therefore, in moving the IF from R to R + dR, enough photons must have arrived at the 
front to ionize all the neutral atoms lying between these positions For unit area of IF, the 
relationship dt= n dR must be satisfied. We get 

dR _ J 

~dt ~ m * 
The velocity given by this expression is the velocity with which the IF moves relative to 
the neutral gas and in the fixed frame of reference (since the gas around the IF is at rest). 

In order to relate J to S* , two effects need to be considered. Firstly, the IF is really 
a spherical surface centered on the Sun, and radiation field at the IF is diluted because 
of spherical geometry. Secondly, recombination takes place continuously inside the ionized 
region and neutral atoms are continuously being created. These atoms absorb photons 
traveling outwards from the Sun and cause a further reduction in the flux at the IF. If an e- 
and ion recombine to form an atom, a new ionizing photon will be required to separate the 
two particles. The photons never reach the boundary R and it will therefore not contribute 
to growth of the region. Ground state recombination is assumed to be balanced by an 
immediate ionization (called the on-the-spot assumption). Hence, 

S = 4nR 2 J + (^j nR 3 n 2 f3 2 (T e ) (2) 

where, f3 2 (T e ), the recombination coefficient into level 2 is given by 

/3 2 {T e ) = 2 10- 6 Tt 3/4) m 3 s- 1 . 
Using equations 1 and 2, the IF velocity is given by 

^ = ^R 2 n l -Rn o p 2 (3) 
dt 4*7T 3 

The IF velocity decreases with increasing R. If the Stromgren radius R s is defined as 
that radius at which the stellar photon output rate just balances the recombination rate in 
the entire ionized volume, then S* and R s are given by 



KR A n^ 2 {T e ) (4) 



35* 2 , 



1/3 



Rs = { ^ 2 p 2 ) (5) 



The characteristic time, tn, for which hydrogen recombines, is 



[n fh ) 
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If the mass surrounding the Hj region is M, the mass per unit area at the separating 
surface will be M/(4ir R 2 ) and the pressure inside will be 2mkT. T is the temperature of 
the ionized region and the factor assumes that the number of ions ni is equal to the number 
of electrons. Neglecting the small gas pressure on the outside of the sphere, we obtain the 
outward acceleration of the boundary as 

(PR _ 2n,kT 

~d~P ~ M/(4irR 2 ) ( ' 

Integrating equation 6 we get a developmental time scale of order 

/3M , ^ 1/2 
= «._•* 



t d = ^—n lk Tj (8) 

To solve Equations 4 and 5, we define the following dimensionless quantities. 

A _ R . Vh _Rs. r _ t dX _ (dR/dt) (g) 
R s tn tji dt Vr 



Equation 3 then becomes 
which has as its solution, 



^ = ^1 (10) 
dt 3A 2 1 ' 



A=(l-e- T ) 1 / 3 (11) 
In the first stage, as shown in fig. 3, R < Rs, Equation 10 or 11 describes the motion. 



3.2 Implications 

We now consider the physical implications of the interactions discussed above. 

1. Using characteristic values only when t=4 wc get A =0.99. Hence the equilibrium 
Stromgren radius is reached only after a long time. 

2. Until R is very close to R s , the radius of the ionized region increases at a speed much 
greater than the characteristic speed (the sound speed, c s w 10 km/s) at which the 
ionized gas reacts to the sudden pressure increase caused by photo ionization. The ion- 
ized gas cannot move appreciably during this phase and the gas density consequently 
is unchanged by the ionization process (hence n^ = n , where n^ is the ionized gas 
density) . 

3. When a time is reached greater than a few times t#, the IF slows down very rapidly 
till its velocity is below the sound speed in the ionized gas. 



3.3 Ionization Heating Equilibrium: 

A shock wave is set up when a piston is advanced supersonically into a gas. The expansion 
velocity of the ionized gas sphere (which is bounded by the IF) is originally equal to about 
Cj, the speed of sound in ionized medium. This is highly supersonic w.r.t the sound speed 
c n (ci/c n < 10) in the neutral gas. The ionized gas sphere plays the role of a piston (fig. 
3) and pushes a shock wave into the neutral gas. The shock then sets the neutral gas into 
motion outwards and it is therefore displaced from its original position. The problem can 
be analyzed using the following assumptions: 
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Figure 3: Interaction of Stellar Wind with ISM 



1. The layer of shocked neutral gas is thin because of the compression across the shock. 
A single radius R and velocity dR/dt then refers to both shock and IF. 

2. The pressure behind the shock wave is assumed uniform in both the neutral and ionized 
gas. 

3. The shock is strong, i.e., it moves highly supersonically into the neutral material. 

4. The recombination rate in the ionized gas balances the stellar UV output rate. 

5. The neutral gas ahead of the shock is at rest. 

Under these conditions, the pressure behind the shock wave, P s , is equal to the pressure 
Pi in the ionized gas given by 

P t = 2n t kT t - n t m H c 2 (12) 

The pressure P s behind a strong shock is related to the shock velocity V s and density 
no by assuming that the shock wave is isothermal. 

Ps = r V s 2 (13) 

Equation 13 can be written as 

P s = Pi=n m H (dR/dt) 2 (14) 
where dR/dt = V s . From equations 12 and 14, neglecting thermal pressure, we get 



dR \ 2 ( rii 



dt I \n 



(15) 



Since the recombination rate in the ionized gas balances the stellar UV output rate 
outside Stromgren radius, using equation 4, 5 and 15, we get 

R 3/2 (dR\ 2 2 / 35* \ V2 = c 2 R 3/2 m 
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Taking the dimensionless variables, 

A = R/R S ;N = c t t/R s (17) 

Equation 16 can be written as 

The boundary condition to solve 18 , is that the time taken to set up the initial Stromgrcn 
sphere is a very small fraction of the lifetime of the H// region, that is N=0 at A=l. Solution 
for equation 18 is 



77V \ 



4/7 



A-(l+— j (19) 
dX ( ' 7N^ ~ 3/7 



Equation 20 implies that dR/dt = Ci at N = 0. 

The condition of final pressure equilibrium is when the hot ionized gas reaches pressure 
equilibrium with the surrounding cool neutral gas, that is 

2n f kTi=n kT n (21) 

Nf =ionized gas density at this stage. Ti and TV are respectively the ionized and neutral 
gas temperatures. The ionized gas sphere must still absorb all the stellar UV photons. Thus 
taking the clue from the equation 4, we get, 

S* = (4/3)irR 3 f n}/3 2 (22) 

Rf =final radius of ionized gas sphere. 
From equations 5, 21 and 22, we get 

(n f /n ) = (T n /2T<); R f = (2T t /T n ) 2 ^R s (23) 

The ratio of fully ionized gas (M/) to the total gas (M s ) contained within the initial 
Stromgren sphere of radius Rs is 

M f /M s = (R f /R s ) 3 (n f /n ) = 2T t /T n (24) 
Typically, n//no ~ 0.005, Rf/R s ~ 34, Mf/M s w 200. Hence at equilibrium we have, 

A = R/Rs w 34 (25) 

Equation 19 then gives corresponding value of N at about 273 and therefore t sa 273 
Rs/cj. 

It is interesting that since t eqm is about 8 10 7 years, typically an order of magnitude 
greater than phase of existence of massive star, it is unlikely that the final pressure equilib- 
rium configuration can be reached unless n is extremely high (n > 3 10 9 m~ 3 ) 
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(d) 



Figure 4 



Regions of flow pattern 



3.4 Interaction of stellar wind with ISM 

The sun emits a lot of gas in the form of a solar wind. The wind velocity (V sw s=s 400 
km/s) pushes the interstellar gas (where the sound speed is d ~ lOkm/s) at a speed which 
is highly supersonic. The wind acts as a piston and a shock wave must immediately be set 
up in the interstellar gas. In doing this, the wind itself must be slowed down. To convert 
the kinetic energy of the wind to its thermal energy, we introduce a second shock into the 
wind itself. There are four distinct regions of flow and three distinct boundaries as marked 
in fig 3.3. 

Region (a) is occupied by un-shocked stellar wind gas moving with a velocity V sw (~ 
400 km/s). It enters the shock SI that converts part of the energy of this wind into thermal 
energy. Region (b) contains shocked stellar wind gas. Since the Mach number of the shock 
is very high, this shock is very strong (see e.g. Harwit, 1998). The immediate post-shock 
gas temperature is 



It is the expansion of this hot bubble of gas, which drives the shell of shocked interstellar 
gas (region(c)), and this is sometimes referred to as an energy driven flow. The expansion 
velocity of the bubble drops because of the loss of internal energy as it does work on the 
surrounding material. This velocity very quickly drops to a value much less than the internal 
sound speed. 

Region (c) is a shell of interstellar gas, which has passed through the outer shock S2. 
Radiative cooling is now extremely effective in this region. Since the shock is a strong 
shock, the upstream to down stream velocity ration (shock strength) is about 4. Using 
the continuity condition, we then get the density immediately behind shock S2 to be four 
times the ambient density. The interstellar gas enters shock S2 at a much lower speed 
than the velocity at which the stellar wind gas enters shock SI. The immediate post- 
shock temperature is therefore considerably lower than that behind SI. The gas behind S2 
contains ions (such as ++ ) which can be collisionally excited by electrons to produce the 
forbidden lines. Cooling is therefore very effective. Region (d) is the surrounding ionized 
interstellar gas with sound speed (10 km/s). Region (c) will be assumed to be a very thin 
shell and for the outer radius of the bubble of hot shocked wind gas. Further, since region 
(b) is thick, we will assume that the shocked wind gas occupies the entire volume interior 
to the shock surface of radius R. Since the interstellar gas is at rest, dR/dt represents both 
the shock velocity relative to the interstellar gas and the expansion velocity of the bubble. 
Physically the thin shell (region(c)) is compressed and accelerated by the hot bubble, which 
has uniform pressure P. 



T — 



32k 



4 10 5 K 



(26) 
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The equation of conservation of momentum of the shell has the form 



d_ 
~dt 



4 d3 dR 
3 at 



= 4nR 2 P (27) 



The left hand side of the equation 27 represents the case of change of momentum of the 
shell, while the right hand side represents the force acting on it. Equation 22 can be written 
in the form 

1 dR 2 ( dR 2 \ 
P = 3* ~t + V°^ R ) (28) 

The thermal energy per unit volume of a mono-atomic gas is equal to 3/2 times the gas 
pressure. The total thermal energy in region (b) therefore equals 4/37ri? 3 (3/2)P. Conser- 
vation of energy for this region is given as: 

dt V 2 J dt dt K 3 ' ' ' 

Equation 29 states that the rate of change of thermal energy of region (b) (figure 3.4is 
equal to the rate at which the energy is fed into the region by gas entering shock SI minus 
the rate at which the hot gas does work on the interstellar gas. This assumes that all the 
kinetic energy of the gas entering shock SI is converted immediately into thermal energy 
Equations 25 and 26 can be combined to give the equation of motion of the shell as 



\ dt 3 J \dt J \dt 2 J V dt 

3 dE* 



2irp dt 



(30) 



We assume the solution to be 



R = At a (31) 



where A and a are positive constants. R then increases with time and tends to zero at 
t=0. Inserting equation 30 in equation 31, we get 

A 5 [a{a - l)(a - 2) + 12a 2 (a - 1) + 15a 3 ]^ 5 "- 3 ) 

= JL(dE*/dt) (32) 
27rp 

Since the mechanical output rate of the star is assumed constant in time both the R.H.S. 
and L.H.S of (29) are time independent and hence a — 3/5. It immediately follows that 

A =fj^y'7^y /5 ( 33) 

Therefore, we find 



1547T J \ Po 



„ ( 125 \ 1/5 fdE*/dt\ 1/5 3/5 

^fi^y^ 1 — y A V 2/s 05) 

dt \154wJ \p dt J V ' 
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Clearly, R and (dR/dt) are related by 

In the equilibrium condition, when the pressure of the solar wind is balanced by the ISM 
pressure, the condition can be specified as 

Psw = Pism (37) 

The SW input pressure is E* . Assuming adiabatic expansion the pressure falls as (see 
e.g. equation 11.31 Longair (1997)). 

^4/3 



E*(r)=E*(0)^—j (38) 
Combining equations 33 and 34 and assuming PISM =1.5 n ism kT ism , at equilibrium, we 

get 

1 .hn sw kT sw {R/Rg) I — \.hrii sm kTi sm (39) 



Therefore, 



Where R„ = Radius of the Sun. 



/?=( ) 3/4 R s (40) 



3.5 Modulation of cosmic rays 

Cosmic rays are charged particles that come to us from the galaxy. The interaction re- 
gion between the ISM and the Sun heavily modulates them. We assume that the generic 
interstellar spectrum of the cosmic rays is given by the function of the form 

f(E) = CE- 2 - 5 (41) 

where E is the kinetic energy of the cosmic rays. This spectral form is assumed based on 
very high energy (> 10 14 eV) cosmic ray data. As we show below, these cosmic rays are 
not significantly modified by the various processes considered here. 

Since these are charged particles, they are significantly affected by the potential difference 
produced due to the interactions discussed above. The potential of interactions is empirically 
defined by 

R 

(eCv sw /k d )dr (42) 

r(E) 

where C is the Compton Getting factor, e the electron charge kd is the diffusion coefficient, 
is given by 

k = (43) 
c 

and 

A = l/(mrV2d 2 ) (44) 

where A is the mean free path, d is the atomic diameter (10~ n m) and n is the density 
of the medium, (see e.g. Lockwood and Webber, 1995, Garcia Munoz 1975, Garcia Munoz 
et al., 1990). The differential flux of cosmic rays at any energy E of cosmic rays is given by 
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Tabic 2: Physical Parameters of different ISM conditions 



Parameter 


Fluff 


Fluff2 


Bubble 


Inter arm 


GMC 


Arm 


n (p/m rf ) 


8.0 10 4 


1.6 10' 


4.0 10 4 


5.0 10 2 


1.0 10 1U 


1.0 10 4 


T e (Kelvin) 


7.0 10 3 


1.0 10 3 


1.0 10 6 


1.0 10 2 


1.0 10 1 


1.0 10 3 


Time (seconds) 


1.0 10 2 


1.0 10 2 


1.0 10 s 


1.0 10 7 


1.0 10 8 


1.010 9 


Recomb coeff. /?2 (m 3 /s) 


2.61 lO" 19 


1.12 10~ 18 


6.32 10~ 21 


6.32 10~ 18 


3.56 10~ 17 


1.1210" 18 


Rccom Fact a 


6.22 10~ 13 


1.64 10~ 12 


5.20 10~ 14 


5.20 10~ 12 


1.64 10" 11 


1.64 10~ 12 


Recomb time t^j (seconds) 


4.78 10 13 


1.11 10 13 


1.98 10 15 


1.98 10 12 


3.51 10 11 


1.11 10 13 


Mism (Kg) 


1.86 10 21 


3.73 10 23 


9.32 10 2() 


1.16 10 19 


2.33 10 26 


2.33 10 20 


V 30und (m/s) 


1.21 10 4 


4.56 10 3 


1.44 10 5 


1.44 10 3 


4.56 10 2 


4.56 10 3 



Tabic 3: Interaction parameters for Stromgrcn radius 



Parameter 


Fluff: 


Fluff2 


Bubble 


Inter arm 


GMC 


Arm 


Stromgrcn radius R s (meters) 


1.13 10 13 


2.02 10 11 


6.18 


10 13 


1.15 10 14 


8.76 10 8 


2.77 


10 13 


Stromgrcn time T s (sec) 


7.38 10 8 


3.51 10 7 


3.39 


10 8 


6.29 10 1() 


1.52 10 6 


4.80 


10 9 


Stromgren vel V s (m/s) 


1.53 10 4 


5.77 10 3 


1.82 


10 5 


1.82 10 3 


5.77 10 2 


5.77 


10 3 


Stromgrcn Temp (K) 


2.10 10 4 


2.10 10 4 


2.10 


10 4 


2.10 10 4 


2.10 10 4 


2.10 


10 4 


Density (/m3) 


2.67 10 4 


7.62 10 5 


1.90 


10 6 


2.38 10° 


4.76 10 6 


4.76 


10 2 


dR/dt (initial) (m/s) 


7.58 10~ 2 


1.21 10° 


5.21 


10~ 3 


1.21 lO" 1 


1.04 10 2 


1.04 


io- 1 


dR/dt (final) (m/s) 


1.12 lO" 15 


5.55 10- 15 


0. 





5.27 lO" 16 


9.81 10~ 13 


0. 





Particle Velocity (m/s) 


2.29 10 4 


2.29 10 4 


2.29 


10 4 


2.29 10 4 


2.29 10 4 


2.29 


10 4 


Mean free path (m) 


8.44 10 16 


2.95 10 15 


1.18 


10 15 


9.45 10 20 


4.73 10 14 


4.73 


10 18 


Diffusion Coefficient (/m 2 s,c=l) 


1.19 10 16 


4.18 10 14 


1.67 


10 14 


1.34 10 20 


6.68 10 13 


6.68 


10 17 


<t> (MV) 


8.64 10 1 


4.44 10 1 


3.39 


10 4 


7.86 10~ 2 


1.20 10° 


3.79 


10° 



where A, B, C, D are constants and equal to 5 108,1150, 650 and 1500 respectively. We 
obtain their values from the current best fits for cosmic ray spectra. E is the kinetic energy 
in MeV/n, E is the rest mass energy, and a is the slope of the incident particles, m=B- 
Ce(-E/^>. 



4 Calculations of the interaction parameters for the Sun 
under varying conditions of the ISM 

There is only one free parameter the Solar UV flux (S*). For the present calculation, we 
take the UV flux to be 10 31 photons/sec (see e.g. Harwit, 1998, Stix, 1989) as the mean 
flux below 90 nm (13.6 eV) to evaluate the ionising strength of the Solar UV radiation. 

We calculate the sizes of the ionization sphere, the ionization pressure equilibrium sphere 
and the kinetic sphere produced around the sun due to different conditions of the interstellar 
medium. We define the various parameter values for the calculations of the interaction 
regions under the different conditions in table 1. We define six specific ISM environments 
as characteristic of regions encountered by the sun during the last 40 million years. Table 
2 gives the values of the constants used and Tables 3, 4, and 5 the derived values of the 
parameters in the interaction regions. We define 2 types of fluff regions (fluff and fluff2). 
The parameters of fluff are taken from the observationally defined values (Frisch, 1998) and 
Fluff2 defines the parameters if the sun were in equilibrium with the fluff. The time scale 
for equilibrium is large and hence we use only the fluff as representing the current situation. 
The difference between the two versions is the mean density and temperature of the LISM. 

The basic features of the calculations are given in figures 5 and 6. The figures show the 
equilibrium interaction radius and the time scale required for reaching the equilibrium for 
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Tabic 4: Interaction parameters for Ionization equilibrium 



Parameter 


Fluff: 


Fluff2 


Bubble 


Inter arm 


CMC 


Arm 


R ther IH (R/H = 34R S ) (m) 


3.83 10 14 


6.88 la 1 '-" 


2.10 lO 15 


3.90 10 i; ' 


2.98 10 1U 


9.41 10 14 


tthermailH (t/ff) (s) 


2.55 10 11 


1.21 10 10 


1.17 10 11 


2.17 10 13 


5.24 10 8 


1.66 10 12 


V/ff (m/s) 


1.50 10 3 


5.68 10 2 


1.80 10 4 


1.80 10 2 


5.68 10 1 


5.68 10 2 


Temperature (K) 


7.00 10 5 


1.00 10 5 


1.00 10 8 


1.00 10 4 


1.00 10 3 


1.00 10 s 


Density (/m 3 ) 


4.00 10 2 


8.00 10 4 


2.00 10 2 


2.50 10° 


5.00 10 7 


5.00 10 1 


dR/dt (m/s) 


6.03 10 1 


2.28 10 1 


7.21 10 2 


7.21 10° 


2.28 10° 


2.28 10 1 


Magnetic field (Tesla) 


1.50 10~ 12 


1.50 10~ 12 


1.50 10~ 12 


1.50 10~ 12 


1.50 10~ 12 


1.50 10~ 12 


Particle Velocity (m/s) 


1.32 10 5 


4.99 10 4 


1.58 10 6 


1.58 10 4 


4.99 10 3 


4.99 10 4 


Mean free path (m) 


5.63 10 18 


2.81 10 16 


1.13 10 19 


9.00 10 20 


4.50 10 13 


4.50 10 19 


Diffusion Coefficient (/m 2 s,c=l) 


7.95 10 17 


3.98 10 15 


1.59 10 18 


1.27 10 20 


6.36 10 12 


6.36 10 18 


<t> (MV) 


2.54 10 2 


3.46 10 2 


8.34 10 3 


1.94 10° 


9.35 10 1 


2.96 10 1 



Table 5: Interaction parameters for Hcliopausc 

Parameter Fluff: Fluff2 Bubble Inter arm CMC Arm 



Rhp (m) 


4.12 10 12 


3.34 10 11 


1.68 10 11 


4.49 10 io 


8.44 10 iu 


8.44 10 i3 


tup (s) 


3.42 10 8 


7.32 10 7 


1.16 10 6 


3.11 10 12 


1.85 10 8 


1.85 10 10 


V H p (m/s) 


1.21 10 4 


4.56 10 3 


1.44 10 5 


1.44 10 3 


4.56 10 2 


4.56 10 3 


Temperature (K) 


6.00 10 6 


6.00 10 6 


6.00 10 6 


6.00 10 6 


6.00 10 6 


6.00 10 6 


Particle Density (/m 3 ) 


1.20 10 5 


3.43 10 6 


8.58 10 6 


1.07 10 1 


2.15 10 7 


2.15 10 3 


Mean free path (m) 


1.87 10 16 


6.56 10 14 


2.62 10 14 


2.10 10 2() 


1.05 10 14 


1.05 10 18 


Diffusion Coefficient (/m 2 s,c=l) 


2.65 10 15 


9.27 10 13 


3.71 10 13 


2.97 10 19 


1.48 10 13 


1.48 10 17 


<t> (MV) 


2.49 10 3 


5.76 10 3 


7.24 10 3 


2.42 10 2 


9.11 10 3 


9.11 10 2 



the Stromgren Radius (SR), the Ionization Front (IF) and the Heliopause (HP). Both the 
axes are in log scales. In the history of the Sun (table 1) indicates that the X-axis indicates 
time before present in the movement of the Sun in the Galaxy. 

The consequences of this changing ISM environment can be quite varied. The present 
studies indicate, for example that when the solar system was passing through a high density 
cloud, the drag on meteors and potential instabilities could trigger a large number of events 
such as increased cometary showers, changes in the solar constant etc. since it is estimated 
that almost 98 % of the interplanetary material is of ISM origin (Frisch, 1998). Also, all 
models of stellar evolution essentially assume that the environmental pressure on the star is 
negligible but clearly, we know that the Heliosphere size is significantly altered by a factor of 
2 or larger even due to small changes in the solar activity (see e.g. Lockwood and Webber, 
1995). The anomalous cosmic rays in the energy range of 10 to 100 MeV/n are particularly 
sensitive to the heliosphcric condition. Since their origin is related to the ionisation of 
interstallar neutrals into the solar system (Cummings and Stone, 1998). Recently, several 
authors: Czechowski et al.,2001; Izmodenov, 2001, Donato, Manrin and Tailler, 2002; and 
references therein. See also Jokipii et al., 1995 and Jokipii, 1996) have attempted to use 
the warm neutral ions as well as galactic cosmic ray modulation to determine the local 
conditions. All these studies rely on measurements of either anomalous cosmic ray studies 
or on galactic cosmic rays in the energy range up to 10 10 eV to determine the current local 
ISM conditions. These are either restricted to deriving inferences of the Heliopause or on 
time average studies of cosmic rays confined in the local bubble. Since the life time of cosmic 
rays is large (15 million years, see Mcwaldt et al., 2001), such studies will not be able to 
give a time sequence of LISM in the past. For high energy cosmic rays however, the origin 
itself is in supernova shock fronts (see e.g. Lagage and Cesarsky, 1981) the gyro radius of 
higher energy cosmic rays may be more sensitive to the LISM conditions since their fluxes 
would be directly modulated by the LISM. 
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Interaction of Sun with ISM 
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Figure 5: Equilibrium interaction radius under different ISM conditions 
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Figure 6: Evolution of the ionization induced interaction regions 



5 Discussion 

We have displayed in figures 5 and 6 the depth of the interaction region under different ISM 
conditions and different physical processes. As can be seen from figure 5 the SR is 4 orders 
of magnitude larger than the HP that is normally considered the boundary region of the 
solar interaction with the interstellar medium. While the SR decreased by about two orders 
of magnitude as the sun entered the local fluff inside the supernova, the HP remained more 
or less constant during the movement of the Sun from bubble to fluff. This is because the 
HP depends on the density and temperature of the medium while SR depends only on the 
ISM density. It should also be noted that the variations in all the three radii can be of the 
order of 6 orders of magnitude. The time scales to reach equilibrium are also quite varied 
(figure 6). As can be seen from the figure, in the low-density region of the bubble, the HP 
reaches its equilibrium in less than 1 year but it would take about 100 years to reach an 
equilibrium in the fluff. On the other hand, the local environment of the GMC, the SR and 



15 



Variation of CR Flux due to Heliopause 
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Figure 7: Variation in the cosmic ray flux due to Heliopause 
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Figure 8: Variation in the cosmic ray flux due to Stromgren Radius 



IF would not reach their equilibrium state for any realistic time scales. 

In order to look for possible direct, observational data to study and understand the 
changing LISM conditions over long time scales, we calculate the modulation of cosmic rays 
over these time scales based on the formulation given above. 

In figures 7, 8 and 9 we have plotted the cosmic ray modulation under different LISM 
conditions. In all cases we assume that the cosmic rays with a spectral index of -2.5 in 
energy (equation 41 impinge on the sphere. We calculate the potential generated by each 
sphere and evaluate the flux parameters. This is the most dramatic of all variations except 
when the solar system enters the bubble. 

In figure 7, the variations due to the changes in the Heliopause alone are shown As can 
be seen from the figure, in the low energy can be as much as 6 orders of magnitude due to 
the variation of the HP alone. Note that a very high energy density in cosmic rays comes 
from the low energy anomalous cosmic rays, which are accelerated in the HP itself. For the 
present study we have ignored this component. 
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Variation of CR Flux due to lonisation Radius 
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Figure 9: Variation in cosmic ray flux due to Ionization front 
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Figure 10: Variation in cosmic ray flux due to variation in the LISM conditions 



In figure 8 we show the effects of SR on the cosmic ray modulations. Since this shock 
front is generated due to ionization and since its radius is much larger, the modulation is 
also more dramatic. Figure 9 shows the modulation produced by the IF which extends 
beyond SR and therefore also affects the CRs more severely. 

In figurelO we have plotted the cosmic ray energy density over a period of time. The 
flux is calculated based on the assumption that the solar system is embedded in a specific 
environment as given by table 1 and that all the three forms of interactions discussed above 
are in equilibrium. We plot the minimum flux that will pass through the three interaction 
regions and take that as the final flux that would be seen by the solar system. Note that the 
two phases of passage through a SN and an Hjj front are not plotted since these are period 
of local and intense cosmic ray re-acceleration which have been ignored in the present work. 
If the present analysis is borne out by further studies, this may be an important component 
in producing the knee in the cosmic ray spectra. 

Florinski Zank and Pogorelov (2003), Florinski, Zank and Axford (2003) and Florinski 
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and Jokipii (2003) have used similar transport calculations to discuss the affects of dense 
molecular cloud on cosmic ray propagation both, for galactic cosmic rays and for anoma- 
lous cosmic rays. Their results show that cosmic rays would be significantly modulated in 
the Heliosphere (Florinski and Jokipii, 2003). The study reported in Florinski, Zank and 
Pogorelov (2003) show work out the modulation formulation for cosmic rays under such 
environment while Florinski, Zank and Axford (2003) show that the 1 0Be data can best be 
explained by assuming that the solar system passed through a molecular cloud 35,000 to 
60,000 years before present. 

6 Conclusions 

We have investigated the interaction of the Sun with the local interstellar medium. We 
show that there have been dramatic changes in the LSIM condition in the last 40 million 
years or longer. We then investigate the physical processes through which the Sun interacts 
with the LISM and show how these regions of interaction would be dramatically different for 
varying ISM conditions. We investigate radiative modulation, radiation induced ionization 
front and the classical solar wind interaction of the ISM. We then investigate the affects 
of these interaction regions on the modulation of cosmic rays. We find that the first two 
parameters will have significant affect on cosmic ray modulation though this has generally 
been neglected in the literature. 
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